Stability of Quasicrystals Composed of Soft Isotropic Particles 
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Quasicrystals whose building blocks are of mesoscopic rather than atomic scale have recently 
been discovered in several soft-matter systems. Contrary to metallurgic quasicrystals whose source 
of stability remains a question of great debate to this day, we argue that the stability of certain 
soft-matter quasicrystals can be directly explained by examining a coarse-grained free energy for a 
system of soft isotropic particles. We show, both theoretically and numerically, that the stability 
can be attributed to the existence of two natural length scales in the pair potential, combined 
with effective three-body interactions arising from entropy. Our newly gained understanding of the 
stability of soft quasicrystals allows us to point at their region of stability in the phase diagram, and 
thereby may help control the self-assembly of quasicrystals and a variety of other desired structures 
in future experimental realizations. 
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Quasicrystals are more common than one had origi- 
nally expected when their discovery vifas first announced.^ 
More than a hundred difi^erent metallic alloys are known 
to form stable quasicrystalline phases of icosahedral 
symmetry alone,'^ with a few dozen additional sta- 
ble phases exhibiting decagonal (10-fold) and possibly 
other symmetries.'^ Yet, to this date, there is no gen- 
eral agreement regarding the origin of their stability 
and the respective roles of enercy and entropy in de- 
termining the observed phases. These growing num- 
bers of stable solid-state quasicrystals, whose build- 
ing blocks are on the atomic scale, have been joined 
in recent years by a host of soft-matter systems ex- 
hibiting qnasiperiodic long-range order with building 
blocks on a much larger scale of tens to hundreds of 
nanometers — micelle-forming dendrimers,- -' star block 
copolymers,'^ mesoporous silica,^ and binary systems of 
nanoparticles.™ These newly discovered soft quasicrystals 
hold the promise for applications based on self-assembled 
nanomaterials,^^ with unique electronic or photonic 
properties that take advantage of their qnasiperiodicityl^-' 
At the same time, they provide alternative experimen- 
tal platforms for the basic study of qnasiperiodic long- 
range order, and offer the opportunity to study the ther- 
modynamic stability of quasicrystals from a fresh view- 
point. To this date, soft quasicrystals have been observed 
only with dodecagonal point-group symmetry, having 
qnasiperiodic order in the 12-fold plane and periodic or- 
der normal to the plane, whereas dodecagonal solid-state 
quasicrystals are rare and mostly only metastable™ Soft 
quasicrystals may belong, therefore, to a distinct class 
of quasicrystals, whose source of stability is likely to be 
different from their solid-state counterparts. We propose 
here a simple theoretical framework to address these new 
systems. We use it to explain the stability of the ob- 
served structures and indicate the (surprisingly simple) 
minimum conditions under which quasicrystals could be 
stabilized. Knowledge of these conditions gives ns the 
ability to estimate the location of the region in the phase 



diagram where quasicrystals should be stable, and thus 
may help control the self-assembly of quasicrystals and 
other desired structures in future experimental realiza- 
tions. 

Several microscopic models have been studied over the 
years, mainly using computer simulations, but also us- 
ing sophisticated analytical methods such as thermody- 
namic perturbation theory,121 to explore the structures 
arising from pair potentials that possess more than one 
microscopic length scale. These studies have yielded 
surprisingly rich phase diagrams even within the lim- 
ited scope of single-component systems, interacting via 
isotropic pair potentials, ^'^"'^ in some cases even find- 
ing stable quasicrystals.^^ is 21 q^^ ^-^^ other hand, phe- 
nomenological models based on coarse-grained free ener- 
gies have been widely applied to treat phase diagrams 
and transitiongSSl and to explain the stability of differ- 
ent phases, including quasicrystals.!^ This is especially 
true in the case of soft-matter systema^SESl j^g ^^ their 
intermediate mesoscopic building blocks, which are sig- 
nificantly larger than the atomic scale, rendering a long- 
wavelength gradient expansion a valid approximation. 

A particular free energy of this sort, which is rele- 
vant for what follows below, was developed by Lifshitz 
and Petrich,'2Sj (henceforth LP) who extended the Swift- 
Hohenberg equatiorp3 to study parametrically-excited 
surface waves (Faraday waves) ^ exhibiting dodecagonal 
qnasiperiodic order. The LP free energy has the form 
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where V^ — d^ + d^ is the 2-dimensional Laplacian. 
It is quite generic, imposing only two requirements on a 
material described by a 2-dimensional density p[x, y): (a) 
the existence of two characteristic length scales, whose ra- 
tio is given by the parameter q; and (b) effective 3-body 



interactions, weighted by the parameter a, that act to 
stabihze structures containing triplets of density modes 
with wave vectors adding up to zero. LP showed that if 
q is chosen around 2cos(7r/12) = v 2 + -x/S ~ 1.93 one 
can obtain a quasiperiodic ground state with dodecagonal 
symmetry, yet no choice of q yields globally-stable ground 
states with octagonal or decagonal symmetry, due to 
insufiiciently-many resonant triplets of modes. Inspired 
by this simple result, we conjectured that the existence 
of two characteristic length scales along with 3-body in- 
teractions may constitute the source of stability of soft 
quasicrystals, all of which (to date) are dodecagonaLl^S 
Here we confirm this conjecture by coarse-graining a mi- 
croscopic partition function for isotropic soft particles 
into an effective free energy. In the limit of small de- 
viations away from the uniform, or liquid, phase this 
coarse-grained free energy can be expanded in a power 
series and mapped onto the simple LP form ([I]) , allowing 
us to gain important insight from the simpler LP model, 
and consequently to explain the stability of the observed 
phases using the full coarse-grained free energy. 

Our starting point is the grand partition function for 
a system of particles with pairwise interactions. 
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where {^n}n=i,...,N are the 2-dimensional positions of 
the centers of N particles, U their pair potential, /3 = 
(k^T)~^ the inverse temperature, and /i the chemical 
potential, which determines the mean particle density, c. 
Using standard methods,'^ one can rewrite the partition 
function in terms of collective coordinates — namely, the 
particle density, c(r) = X]n=i <5(r — r„), and its conjugate 
field — rather than discrete positions. At the mean-field 
level, which amounts to a saddle-point approximation for 
the integration over the conjugate field, the transformed 
partition function becomes Z — j 'Dce~^^^'^' , where 



■^ [c(i-)] = ^ /" drdr'c(r)C/(r - r')c(r') 



+ / dr{fcBrc(r)[lnc(r) - 1] - /ic(r)}. (3) 

The coarse-grained free energy functional given in Eq. ([3]) 
contains the familiar mean-field terms of pair interaction 
and ideal entropy. Although it could have been written 
from the outset, we wish to highlight the ability to extend 
the current theory to higher order, particularly in light 
of the cautionary remarks of Schwartz and Vinograd.!^ 
We shall assume that the equilibrium density field is the 
one that minimizes J-' for the given T and /i and the spe- 
cific choice of C/(r). Results of such direct minimization 
will be presented shortly. However, since it is not a priori 
obvious what pair potentials and thermodynamic param- 
eters may yield quasicrystalline order, it is beneficial first 
to characterize J^ and relate it to J^p of Eq. (fTl) . 



Above a certain critical temperature, T > Tc, the equi- 
librium density for any /x should be uniform, c(r) = c. 
Minimizing T with respect to such a uniform field yields 
the relation. 



/i = kBTlnc + Uqc, 



(4) 



where Uq = J drU{r). For T < T^ the equilibrium den- 
sity is expected to become nonuniform at a certain value 
of fi (or, alternatively, above a certain mean density c). 
Assuming that T is only slightly smaller than Tc, we sub- 
stitute c(r) = c[l -I- p(r)] in J^ and expand to 4th order 
in small p. The result can be written as 
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where tildes denote Fourier-transformed quantities. The 
critical temperature, below which perturbations start to 
grow, is given by 
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where C/min is the minimum of the Fourier-transformed 
pair potential, which must be negative for T^ to be pos- 
itive. The approximate coarse-grained free energy ([5]) 
resembles J^p ^, with the gradient term replaced by a 
Fourier-space integral of the pair potential, and the ex- 
pansion in powers of the density arising from the entropy 
term of Eq. ^. 

To try to predict the conditions under which dodecago- 
nal quasicrystals minimize J^ of Eq. (|3| , let us recall the 
findings of LP, who showed that such crystals are sta- 
bilized in their simpler model, J^lp of Eq. (Il]), under 
the following conditions: (i) q ~ 2cos(7r/12); (ii) < 
e/a^ < 0.088. LP further showed that, for e/a^ > 0.088, 
a hexagonal state is obtained. In terms of our free en- 
ergy expansion (l5| , condition (i) translates into a require- 
ment for the properties of f7(k). The first two minima of 
this function should be located at wavevectors ki and fc2 , 
whose ratio q = ^2/^1 is around 2cos(7r/12) ~ 1.93, and 
should have a similar depth. Condition (ii) translates 
into a requirement for the thermodynamic parameters, 
T and n (or c), namely, < f - T/T^ < 0.066, i.e. T is 
restricted to a small range below Tc- 

We test these estimates, obtained from the approxi- 
mate free energy expansion ([5]), against the numerical 
minimization of J^ [c(r)] of Eq. (Isl). We find that a vari- 
ety of isotropic pair potentials, three examples of which 
are shown in Fig. Ill contain a sufficient number of tun- 
able parameters to satisfy condition (i) above. The form 
of these potentials is motivated by the qualitative fea- 
tures of the experimental systeiUjISI consisting of spher- 
ical micelles whose interaction should contain an inner 
repulsive core, a region governed by van der Waals at- 
traction, and a longer-range soft repulsion. Using these 
potentials in Eq. (|3|, we find the minimum free-energy 
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FIG. 1. Three of the isotropic pair-potentials used in this study (top row), whose Fourier transforms (bottom row) are designed 
so that the ratio q of the positions of their first two minima is approximately 1.93. The first potential in (a), motivated by the 
structure of dendritic micelles, consists of a normalized repulsive core, with a short range (van der Waals) attraction at the 
core boundary, and an exponentially-decaying soft shoulder. This behavior is approximated by even simpler 3-step and 2-step 
potentials, shown in (c) and (e). Potential parameters: (a) U{r) = 1 for < r < 1; and —u\/{r- — 1 + 5) + U2 exp[— (r — 1)/A], 
with Ml = 0.02057, u-z = 1.057, A = 0.9756, and ,5 = 0.01, for r > 1. (c) U{r) = 1 for < r < 1; -0.0770 for 1 < r < 1.0485; 
0.4690 for 1.0485 < r < 1.6570; and for r > 1.6570. (e) U{r) = 1 for < r < 1; 0.457 for 1< r < 1.7442; and for r > 1.7442. 



state by numerically integrating the corresponding relax- 
ational equation, dtC = —6T/Sc, using a pseudospectral 
method — with local terms evaluated in direct space and 
non-local terms in Fourier space — starting with random 
initial conditions and waiting until a steady state is ob- 
tained. We note that this equation does not describe the 
actual dynamics and is used merely as a minimization 
tool. 

When q is selected according to condition (i) — as 
demonstrated in Fig. IT] — and T and fi are set such that T 
is just below T^ [condition (ii)], we indeed find a minimum 
free-energy state which is a dodecagonal quasicrystal, as 
shown on the left-hand column of Fig. [2] fo r the partic- 
ular case of the 3-step potential of Fig. [l|c). When T 
is decreased slightly further a transition is observed to 
a hexagonal state [middle column of Fig. [2], restricting 
the stability of dodecagonal quasicrystals to a narrower 
range than that predicted by the approximate free en- 
ergy expansion ([5]).'^^ The same stabilizing mechanism of 
two length scales and 3-body interactions can be used to 
obtain other structures pS for example, setting q — \/3 
yields immediately below Tc a hexagonal crystal, shown 
on the right-hand column of Fig. [2] We have observed 
that slight variations in the potential parameters still 
yield the expected phases. A systematic study of the 
stability boundaries in parameter space is left for future 
work, perhaps using more realistic potentials, derived 
from specific experimental realizations. 

While the first two minima of the Fourier-transformed 
pair potential U{k) must be negative to give a positive 



Tc, the value of U{k — 0) need not be negative. This 
allows a purely repulsive pair potential U{r) to stabilize 
a quasicrystal, as we show here for the simple two-step 
potential of Fig. W[e). Despite the appealing simplicity 
of this potential — defined after scaling by two parameters 
only — we emphasize that potentials of this sort can sat- 
isfy the requirements on the Fourier space minima only 
near a single choice of their two parameters. As a con- 
sequence, attempting a naive numerical search for a do- 
decagonal quasicrystal with this potential, by setting the 
extent of the second step to i? = 2cos(7r/12) ~ 1.93, 
would fail. We now understand that in order for non- 
linear mode interactions to stabilize a dodecagonal qua- 
sicrystal, one must adjust this ratio not in the real-space 
pair potential, but rather in its Fourier transform. This 
is obtained in this case by setting R ~ 1.74 in real space. 
We have managed to stabilize a dodecagonal quasicrys- 
tal with the simple two-step potential by eliminating the 
need to search around in parameter space. Our theoret- 
ical understanding of the source of stability allows us to 
point at the stability regions in parameter space, even 
when these regions are extremely narrow. 

Thus, we confirm that the existence of two char- 
acteristic length scales and sufficiently strong 3-body 
or higher-order nonlinear interactions can account for 
the stability of dodecagonal quasicrystals of isotropic 
soft particles. More specifically, we determine the two 
length scales through requirements on the minima of the 
Fourier-transformed pair potential and, alth oug h 3-body 
terms may arise from various interactions,'^^ we show 



FIG. 2. Real space densities c(r) (top row) and their Fourier transforms |c(k)|, with c(k — 0) removed (bottom row), obtained 
by a numerical minimization of the free energy ([3|, with the 3-step potential of Fig. [I] Other potentials yield similar results. 
Dodecagonal crystal (left column): 3-step potential, with parameters specified in Fig. [l[c). Thermodynamic variables are 
fi — 2.2895 and T — 0.999Tc. Hexagonal crystal (middle column): Same as in the left column, but with a lower temperature 
T — 0.980rc. Hexagonal crystal (right column): 3-step potential similar to Fig. fll^c), with parameters U{r) = 1 for < r < 1; 
—0.115 for 1 < r < 1.053; 0.370 for 1.053 < r < 1.600; and for r > 1.600, designed so that q ~ \/3. Thermodynamic variables 
are /i = 1.9041 and T = 0.999rc. 



that translational entropy suffices to provide the required 
term. Thus, the delicate interplay between interaction 
and entropy can give rise to stable quasicrystals even 
for relatively simple isotropic potentials. By designing 
these potentials along the guidelines provided here, one 
should be able to control the self-assembly of quite com- 
plex structures. Preliminary tests, using molecular dy- 
namics simulations, indicate that our design principles 
seem to work.'^^ 

This work can be extended in several directions. The 
results of the approximate theory presented here should 
be verified using direct computer simulations. Impor- 
tantly, such simulations can be used to clarify the actual 



stability regions of the quasicrystalline phase. The dy- 
namics of crystallization, and of collective degrees of free- 
dom in the ordered state, can be studied by replacing the 
equation used here to minimize the free energy by one 
that is adequate for a conserved density field, possibly 
while taking thermal fluctuations into account. Finally, 
a similar coarse-graining procedure could be applied to 
two-component systems or anisotropic potentials.'^ 
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